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Abstract 

The Ap interaction close to the £JV threshold is considered. Specifi- 
cally, the pronounced structure seen in production reactions like K~d — s> 
ir~ Ap and pp —¥ K + Ap around the EA threshold is analyzed. Modern 
interaction models of the coupled AN — "EN systems generate such a struc- 
ture either due to the presence of a (deuteron-like) unstable bound state or 
of an inelastic virtual state. A determination of the position of the promi- 
nent peak as observed in various experiments for the two aforementioned 
reactions leads to values that agree quite well with each other. Further- 
more, the deduced mean value of 2128.7 ± 0.3 MeV for the peak position 
coincides practically with the threshold energy of the E + n channel. This 
supports the interpretation of the structure as a genuine cusp, signaling 
an inelastic virtual state in the 3 5i — 3 D\ partial wave of the EAT isospin 
1/2 channel. There is also evidence for a second peak (or shoulder) in the 
data sets considered which appears at roughly 10-15 MeV above the EA 
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threshold. However, its concrete position varies significantly from data 
set to data set and, thus, a theoretical interpretation is difficult. 

1 Introduction 

Since the discovery of strangeness, the hyperon-nucleon (YN) interaction has 
been of fundamental interest both theoretically pQ as well as experimentally [2] . 
First, its knowledge is important for the general understanding of the struc- 
ture of hadrons and their constituents. Further, it is also needed to explain 
the spectra of hypernuclei, which, in return, also convey information on these 
interactions especially at small relative energies. Unfortunately, since hyperons 
are short lived, this energy region is practically inaccessible by hyperon beams. 
Thus, possible bound states are indispensable as the source of information. In 
this report we will concentrate on the Ap interaction at energies close to the 
threshold of the SiV channels. These are at 2128.94 MeV (for S+n) and at 
2130.9 MeV (for £°p). In this region experimental data for elastic scattering, 

Ap -> Ap , (1) 

indicate an enhancement in the Ap cross section as we will discuss in the next 
sections. However, the dynamical origin of this enhancement remains unclear 
so far. It could be a cusp structure due to (and at) the opening of the UN 
threshold and then would be a signal for an inelastic virtual state (we follow 
here the nomenclature used and explained in Ref. [5J) or due to a bound £°p 
or Ti + n state, i. e. a deuteron-like but unstable bound state. In the latter case 
the peak of the cross section has to be below the EAT threshold. In principle, it 
could also be a Ap resonance above the S7V threshold. 

Data from elastic scattering experiments that cover this range exist in the 
literature, but so far the momentum resolution of the A beams has been insuf- 
ficient to draw firm conclusions. A much more promising avenue is offered by 
the study of final state interactions (FSI). Assuming relative weakness of the 
pion interaction one possibility is the strangeness exchange reaction 

K^d^n-Ap. (2) 

Also strangeness production processes like 

tt + d — > K + Ap (3) 
pp -> K + Ap (4) 

should contain basically the same information. 

In this work we will concentrate on the two reactions ([2j and (j4j . With re- 
gard to the former reaction, evidence for an enhancement in the Ap cross section 
near the EAT threshold has been already found in the late 1960s and confirmed 
in later experiments [4TllO|. We review those data and we also re-analyze them 
with the aim to determine accurately the position of this enhancement. Exper- 
imental information about the reaction |4j in the region of the T,N threshold, 
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has become available much more recently and clear evidence for the presence of 
an enhancement at that threshold is only emerging right now. Here we perform 
an analysis of data from the inclusive measurements of the reaction pp — > K + X 
performed at Saclay [TT] and in Jiilich [T2] , respectively, and attempt to extract 
the enhancement and its position from those experiments too. This sort of anal- 
ysis suffers from lacking precise pp exclusive data to rely on. However, as we 
will see, the positions determined from the two reactions ([2]) and (j4]) agree with 
each other, and they also coincide with the opening of the S + n channel within 
the error bars of our analysis. 

We also study an additional peak (or shoulder) that is present in basically 
all the aforementioned measurements and located a few MeV above the EiV 
threshold. However, in this case it turns out that there are sizable variations 
of its position between reaction ([2]) and Q , but even between measurements of 
one and the same reaction. Thus, the physical significance of that peak remains 
unclear to some extent. 

The paper is organized as follows. In the next section we take a look at the 
status of the results from elastic Ap scattering. We discuss also the behavior 
of the Ap cross section around the SiV threshold as predicted by various YN 
interaction models from the literature. In the two subsequent sections the data 
for the two reactions ^ and @ are discussed. In section 5 the peak structures 
found are analyzed. In particular, we determine their position as seen in the 
various measurements for the two reactions in question. The paper ends with a 
summary. 

2 Elastic Ap Scattering 

Before discussing the situation for AN scattering let us briefly recall some well- 
known features of coupled-channels dynamics [3l[T3] . Conservation of flux and 
the associated unitarity of the S'-matrix necessarily imply anomalies at the open- 
ing of new thresholds [13] . Specifically, at an S'-wave threshold the cross section 
of the "old" channel as a function of the energy will, in general, have infinite 
slopes at such a threshold. The resulting structures are usually called cusps or 
rounded steps, depending on their specific shape [3j[l3]. Whether these struc- 
tures remain primarily of academic interest or manifest themselves via large, 
experimentally observable effects depends strongly on the strengths of the inter- 
actions in the coupled channels. In particular, pronounced threshold phenomena 
always go along with near-by poles in the scattering amplitudes of the involved 
channels that are associated with (inelastic) virtual states or (unstable) bound 
states [SHU]. 

Modern meson-exchange models of the YN interaction such as the Jiilich 
[H1[T6] or Nijmegen potentials [El [18] are derived under the assumption of 
(broken) SU(3) symmetry. This symmetry implies, that the strongly attractive 
forces that yield the deuteron bound state (in the 3 Si — 3 Z?i partial wave) and 
a virtual state in the 1 So partial wave in case of the NN system will likewise 
act in the strangeness S = — 1 sector (see, e.g., Refs. [19][20] f° r details on the 
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Figure 1: Elastic A A" cross section as function of the center-of-mass energy. 
Jiil04 and Jiil89 are results for the Jiilich potentials published in Refs. [T5] 
and [16], respectively. NSC97f and NSC89 refer to results of the corresponding 
Nijmegen soft-core potentials [17] and [18]. Results obtained at leading-order 
chiral EFT p~9j[2n] are indicated by the grey band. The dashed line is the 
threshold for the Ap — >• EA" transition. The left part has an expanded energy 
scale. 



SU(3) relations). Specifically, there is a strong coupling between the AN and 
EA" systems. It is caused by the long-ranged tensor force provided by pion 
exchange and boosted by the fact that the thresholds of the two channels are 
only separated by 77 MeV. Therefore, it is not surprising that practically all 
YN interactions that fit the data and include explicitly the coupling between 
the AN and EAT channels predict also sizeable threshold effects. 

A closer inspection of the results for published interaction models reveals 
that, in general, they can be grouped into two categories. In one case the 
predicted AN cross section around the EA" threshold is in the order of 40 to 50 
mb and, more characteristic, roughly a factor four larger than a few MeV away 
from the threshold. The cross section is so large because one of the eigenphases 
of the tensor-coupled 3 5i — 3 D\ partial wave (in most cases the 3 Di) passes 
through 90° right below the EiV threshold. Due to the latter aspect, the peak of 
the cross section is actually not at the EAT threshold but slightly below. Thus, 
in this case the AN cross section exhibits a typical resonance-like behavior. 
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Moreover, no cusp appears at the actual EiV threshold, only a rounded step 
[3lH3]. However, since the peak of the cross section occurs so close to the EA 
threshold - often the separation is less than an MeV - it is usually impossible 
to recognize the above features in the published results due to the scale used 
for the figures! Pole searches performed for the amplitudes produced by the 
potential models in question found that these are located in the second quadrant 
of the complex plane of the relative momentum in the EA channel [14] . Thus, 
these YN interactions are characterized by the presence of an unstable bound 
state, i.e. a deuteron-like EA state [3J. The Nijmegen potentials NSC97f [T7] 
and NF [2T] but also the original Jiilich YN interaction [TH] belong to this 
category. In the case of an unstable bound state it is also possible that the 
3 Si eigenphase passes through 90° (instead of the 3 D 1 ). Interestingly, such a 
scenario is seldom realized. In fact, we are aware of only one meson-exchange 
YN potential where this happens, namely the Nijmegen ESC04 interaction [22] . 
In case of the interactions considered in [23] it is also the 3 Si that passes through 
90° . But since these potentials were intended for application in Faddeev-type 
calculations for simplicity reasons only S-waves were taken into account. One 
of the interactions considered in [23] has the rather unique feature that the 
predicted 3 Si phase passes through 90° slightly above the EA threshold. As 
far as we can see, this does not happen for any of the meson-exchange YN 
potentials whose phase shifts are documented in the literature. 

The second category of YN potentials produces a peak in the AA cross 
section precisely at the EA threshold. Thus, now we do observe a genuine 
threshold cusp. Here none of the relevant eigenphases passes through 90° . In 
general the AN cross section at the EA threshold is roughly a factor two larger 
than at a few MeV away from the threshold. Usually, the peak values are 
around 20 mb, but can still reach up to 40 mb. The poles for this kind of 
potentials are located in the third quadrant of the complex plane of the relative 
momentum in the EiV channel [14] . They are an indication for the presence of 
inelastic virtual states [3J, i.e. the analog of the virtual state in the NN 1 Sq 
partial wave. The Nijmegen potentials ND [23], NSC89 [18. and ESC08 [25], and 
also the recent Jiilich YN interaction [TS] belong to this category. Also a YN 
interaction derived at leading-order in chiral effective field theory (EFT) [H?1I2U] 
predicts such a behavior. 

In order to illustrate the statements made above, Ap cross sections for the 
Nijmegen soft-core potentials NSC97f pj] and NSC89 [18], and for the Jiilich 
potentials from 2005 [15] and 1989 [16] are shown in Fig. [1] for energies around 
the EA threshold for two different energy scales. Results obtained at leading- 
order chiral EFT [19,20 are indicated by the grey band. One can see that the 
cross sections predicted by those YN potentials indeed exhibit different features 
at the EA threshold. As said above the Nijmegen NSC97f and the Jiilich 1989 
one-boson exchange models produce a deuteron-like unstable bound state in 
the EA channel. In both cases the 3 D\ AN phase shift shows a resonance-like 
behavior and crosses 90° slightly below the EA threshold 15, 19 . The rounded 
step in the cross section of the Jul89 potential is clearly visible in the left figure 
(with magnified scale). The structure produced by the NSC97f potential is 
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similar. However, since in this case (in our calculation) the 3 Di phase shift 
crosses 90° at a mere 20 keV below the nominal EiV threshold a further increase 
in the scale would be required to see that. The Jul04 and the NSC89 models 
and the EFT interaction support an inelastic virtual state rather than a bound 
state and, consequently, a genuine cusp structure appears in the cross section. 
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Figure 2: Cross sections for elastic Ap scattering. Data are from [21 I26H30] ■ The 
curves are results from YN models, namely from the Nijmegen YN soft-core 
potentials NSC97f p7] (dashed curve) and NSC89 [18] (dash-dotted curve) , and 
from the Jiilich one-boson-exchange models [15] (solid curve) . and [16] (dashed 
curve). Results obtained at leading-order chiral EFT 19, 20] are indicated by 
the grey band. The thresholds for the reactions Ap — s- TiN are indicated by 
arrows. 



Note that that the results in Fig. [T] are calculated in isospin basis using 
isospin-averaged S and nucleon masses. In particular, the AN momentum and 
the HAT threshold are evaluated for = [m p + m n ]/2 and = [m s + + 
m.£0 + m£-]/3, respectively. 

Experimentally, in principle, the Ap interaction could be studied in the range 
of interest, i. e. in the vicinity of the UN thresholds by elastic scattering. In 
Fig. [2] all cross sections [2 , 26 30 - to the best of our knowledge - for elastic 
scattering in a momentum range from 400 MeV/c to 1 GeV/c are collected. 
In addition to data we include the model calculations for the Nijmegen soft- 
core potentials NSC97f [17 and NSC89 [TS], the Jiilich meson-exchange models 
Jiil05 [15] and Jiil89 [16], and the leading-order chiral EFT interaction [191120] . 

Here the computation of the cross section was done in particle basis so that 
the E + n and E°p thresholds could be correctly implemented. Partial waves up 
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to L < 2 have been taken into account. Note that the agreement between data 
and calculations at low energies (not shown here) is of similar quality for all 
models. 

Obviously, the data on elastic scattering are insufficient in quality to allow 
to discriminate between the different scenarios. It is therefore highly desirable 
to have additional data of high quality. 

3 The reaction K~d — > ir~Ap. 

One possibility to study the Ap interaction is via the FSI in reactions like 
K~d —> ir~Ap. In elastic scattering the contribution from the spin-triplet waves 
to the cross section can be expected to be the larger one due to its statistical 
weight. For the above reaction we expect likewise the main contribution to come 
from spin-triplet states and, specifically, from the 3 Si-wave, but due to different 
reason: the deuteron is already present in the initial state. For kaon absorption 
at rest the final baryonic state has to have predominantly the quantum numbers 
of the deuteron. At lowest order the reaction should be dominated by K~n — > 
7r~A (quasi-elastic) scattering leaving the spin-space part of the baryon state 
unchanged. 

Fig. [3] shows the Ap spectra in the center-of-mass (cm.) system. The data 
are taken from Refs. j3H7,9,10,31J. There is clearly a peak around the UN mass, 
i.e. around 2.13 GeV. However, the shape of the peaks as well as the underlying 
cross section vary from one experiment to the other. One reason for this are 
different boundary conditions in the experiments and the analysis. We will 
return to this point later. In order to study this dependence further we compare 
in more detail some of the spectra. The data from Sims et al. [7] are ignored 
since they have a poor signal to background ratio, i.e. do not allow to study the 
structure in detail. Similarly, the data from Alexander et al. [31] as well as those 
from Cline et al. [5] with bin widths of 10 MeV and poor statistics are excluded. 
The latter data have only two bins with 67 counts together in the peak region 
thus a width of the peak cannot be extracted. The data from Eastwood et al. [9 
on the other hand contain some more counts (130) but in thirteen bins. We 
therefore kept these data. First we inspect the spectra from Tan [4] and Braun 
et al. [6j. These are data sets with rather large statistics. The smooth yield 
below the peaks contains, in addition to the direct reaction contributions 
from the reaction with a heavier intermediate hyperon: K~d — > £~(1385)p 
which decays in a second step as E~(1385) — > Att~ . In order to avoid any 
phenomenological modelling and possible ambiguities associated with that we 
simply subtract some yield below the peak by fitting polynomials to these yields. 
We will discuss this further in section [531 The result is depicted in Fig. 2J Then 
this fitted cross section is subtracted from the experimental cross section so that 
only the peak structure remains. The peak in Tan's paper [3] shows a shoulder 
on the heavy mass side. In the case of the unconstrained data from Braun 
et al. [6] the structure looks quite different. There is much more yield to the 
high-mass side than in Tan's data. At this point we have to elucidate that 
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Figure 3: Invariant mass spectra for the reaction K~d —> ir~Ap. The data are 
from Tan [4] , Cline et al. [5] , Braun et al. [6] , Sims et al. [7] , Alexander et al. [8] , 
Eastwood et al. [9] , and Pigot et al. [10] . The dotted line indicates the averaged 
EiV mass 2.13 GeV. 
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Figure 4: The missing mass spectra for the constrained data (see text) from [6 
(upper panel) and [4] (lower panel). The solid curves show fits to the data where 
the peak region has been excluded. 
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this comparison is made between different things. The data from Braun et al. 
have no constraints while in case of Tan, counts with proton momenta below 
75 MeV/c have been cut. Thus, reactions with the proton as a mere spectator 
K~d — > 7r~ Ap s , with the two baryons being uncorrelated are excluded. In order 
to have really both hadrons in the entrance channel participating in the reaction 
Braun et al. introduced two cuts: (i) a threshold in the proton momentum of 
150 MeV/c and (ii) the requirement of the angle between the incoming kaon 
and the outgoing pion 0.9 < cos(K, n) < 1. If one takes that into account the 
peak agrees, with respect to its shape and its position, to a large extent to the 
one reported by Tan [3] as is discussed in Sect. [S] Note that Eastwood et al. [9 
required the proton momenta to be larger than 170 MeV/c. 

While almost all groups relied on bubble chambers, Pigot et al. [10] applied 
a magnetic spectrometer to detect the emerging pions. However, most 7r~'s 
are not from reaction @ but from beam decays K~ — > 2tt~tt + . In order to 
reduce the number of such events the target cylinder was surrounded by twelve 
scintillation counters. After the decay of the A into two charged particles one has 
then three charged particles in addition to the forward going pion. Therefore, 
the authors studied the data with charged particle multiplicity m > 2 in these 
scintillators. For the width of the peak they give only upper limits. They 
also studied the line reversed reaction ([3]). The cross section for this reaction is 
smaller than for the strangeness exchange. Spectator protons are rarely detected 
in the set-up since they are stopped in the liquid deuterium target being 4 cm 
in diameter. Their final result is 2129±0.2±0.2 MeV and 16.7±1.9±2 MeV 
for the position and the width of the peak, respectively. Here we make use of 
their data set of the reaction ([2]) taken at a beam momentum of 1.4 GeV/c and 
multiplicity m > 3. For this case the spectrum is rather clean and the statistics 
is still sufficient - which is not the case for the other data sets. We then proceed 
as in the other cases. The final results for the structure at the UN threshold 
are given for all cases in Figs. [9land[TT] 

4 The reaction pp — > K + Ap 

Most studies of the reaction pp K + Ap report only total cross sections. Hogan 
et. al [32] and more recently the COSY-TOF collaboration [33,34 published 
spectra. However, these experiments although having rather large acceptances 
suffer from insufficient resolution to study a peak in the threshold region. These 
unfavorable boundary conditions were overcome with sufficient high resolution 
by Siebert et al. [TT] employing the SPES4 spectrometer at Saclay and, more 
recently, even more by the HIRES experiment [12 , 35 making use of the Big Karl 
spectrometer [36] at the COoler SYnchrotron COSY Jiilich. The disadvantage 
of the experiments from Refs. [ITJ[l2l|32]|35] is that they are inclusive. This 
means that above the SAT threshold also the £°p and E + n channels contribute 
to the experimental cross section. Thus a peak in the experimental spectrum 
is more difficult to see because the signal will be distorted by the rising contri- 
butions from the UN channels. Therefore, the following analysis is incomplete 
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and should be revisited upon future availability of good near-threshold E — N 
production cross sections. 

Already at first glance there are differences between reactions © and ((3]). 
For example, available data suggest that the peak to background ratio of the 
structure around the T,N threshold is significantly smaller for the latter reaction. 
Furthermore, the detection of a kaon in the final state is much more delicate 
as compared to the detection of a pion, due to its short life time. One thing 
in common, after applying the cuts as discussed in the previous section, is that 
both baryons participate in the reactions. 

Here we will first concentrate on the data from Refs. QT] and [T3]. In these 
experiments the four momentum of the kaon was measured and the missing 
mass of the YN system was deduced. The missing mass resolution in the Saclay 
experiment [IT] varied between 3 to 5 MeV (FWHM) depending on the angle 
while the HIRES experiments [35] and [T2] reported a resolution of 0.8 MeV. 

The bubble chamber measurements and the exclusive COSY-TOF data are 
projections from Dalitz plots and hence they are in the cm. system. To make 
the other data comparable to them we transform them into the cm. system. 

The authors of Ref. [11] report about several energy shifts they have applied 
to their data. We shift the data back by two MeV in order that the steep 
rise of the cross section agrees with the threshold for the Ap channel. This is 
demonstrated in Fig. [5] where the shifted data are compared with calculated 
cross section 

j^Er * $3 ' Ffsi (5) 

with $3 the three body phase space and Fpsi is the enhancement factor for 
the final state interaction from Ref. [35]. In order to show the effect of 2 MeV 
we have shifted this curve by two MeV up. This curve cannot account for the 
data points in the steep rise close to threshold (see insert in Fig. [5]). The data 
close to the threshold agree well with FSI modified phase space. This Ap yield 
represented by Eq. ([5]) is then subtracted from the experimental spectrum. 
The remaining yield is then assumed to be due to the contributions from the 
pp — > K + Yi + n and pp — > K°T, + p reactions, and from the structure of interest. 
Finally we subtract the cross section of the two E channels assuming K + T,N 
being phase space distributed. This is also shown in Fig. [5J 

The same procedure is applied to the HIRES data. The different steps are 
depicted in Fig. [5j By this method the structure close to the EiV thresholds is 
extracted. 

The new data from TOF [37] are exclusive. In general such data are superior 
to inclusive data since subtraction of cross sections due to other reactions is 
unnecessary. The spectra as function of the Ap mass were obtained by projection 
of the Dalitz plot onto this axis. The particles detected in the final state were 
K + , p and the decay A — > ir~p. The first set of data labelled here by TOF I were 
taken at a beam momentum of 3.003 GeV/c while the second one was taken 
at 2.95 GeV/c. The resolution of the first set is 6 MeV after a kinematical fit 
while for the second 2 MeV is reported due to a new tracking system [38]. We 
proceed as before by fitting Eq. © to the data outside the peak region. This 
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Figure 5: Upper curve: The data from SPES 4 (full dots) [IT] at a lab. angle of 
(10±2) °. The beam had a momentum of 3.1 GeV/c. The data are transformed 
into the cm. system and shifted by 5 nb/MeV sr in height and by 2 MeV down 
in the mass scale. The solid curve is the fit with Eq. (|5|) making use of the 
FSI parameters from Ref. [35] . The lower curve (full triangles) is the difference 
between the data in the upper curve and the fit. The short dashed curve is the 
fit of phase space for S + n. The insert shows again the solid curve (without the 
5 nb/MeV sr shift and the data without the mass shift of 2 MeV.) 



is shown for the TOF II data set in Fig. [7] This procedure accounts very well 
for the exclusive data. Especially the high mass side is well reproduced by pure 
phase space. We can still go a step further and compare the peak regions of 
the exclusive data with the inclusive data. This comparison is made in Fig. |5J 
There is a remarkable degree of similarity between the two results although the 
latter data are taken at a beam momentum 80 MeV/c higher than in the case 
of the HIRES data. The similar shape in the right side of the peak - except 
for m(Ap) s» 2.1475 GeV - contradicts the statement made in Ref. [39 that a 
large fraction of the cross section ascribed to the T,N production might be Ap 
production. 

At this point it might be appropriate to comment on that work further. 
First the authors state that "it is further hypothesised (in Ref. [12]) that, away 
from the near-threshold region, the differential cross section for pp — > K + pA is 
essentially the same to the left and right of the UN threshold, in marked contrast 
even to the phase-space behaviour." Inspection of the figure in Ref. [3H| shows 
that FSI was ignored when showing phase space dependence. This is in stark 
contrast to the experimental finding, as can be seen in Fig. [7J that including 
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Figure 6: Same as Fig. [5] but for the HIRES data [12] measured around zero 
degree and at a beam momentum of 2.87 GeV/c. 



FSI makes the shape of the cross section underlying the peak rather flat. This 
is even more pronounced for the study in the lab. system. 

5 Analysis of the peak structure 

Finally, Breit-Wigner distributions 

d 2 a _ l_da T/2 

dD,K^dm^p 7r df^A'.Tr {m(Ap) — mo) 2 — (r/2) 2 

were fitted to the different yields. Fit parameters were the width T and mo, the 
centroid of the peak. The production cross section da / <KIk,-k is also fitted. In 
Fig. [9] these fits were compared to the data. The deduced centroids and widths 
are shown in Fig. [9] 

The fits yield a mean for the centroids mo = 2130.3 ± 0.6 MeV and for the 
widths r = 12.2 ± 1.3 MeV. There is a large scattering among the values and 
this is reflected by large x 2 /dof values. Especially the results for the data set 
TOF I are larger than the other values. This data set shows an asymmetric 
peak with a centroid at 2.132 ± 0.006 GeV [37] ■ When a symmetric peak is 
fitted the centroid moves further up to 2.1376 GeV. In particular, this is also in 
disagreement with the values extracted for the recent and more precise TOF II 
data. We, therefore, exclude these data from further analysis. 

From inspection of the fits one gets the impression that the peak centroids 
are, with the exception of the data from Braun et al. [6], above the S + n thresh- 
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Figure 7: The exclusive cross sections for the reaction pp — > K + Ap at 2.95 
GeV/c from Ref. [37l[38]. The curves indicate the fit to the data outside the 
peak region. The solid curve is the fit with Eq. ([5]), the dashed curve is the 
pure phase space part. 



old which is also indicated in Fig. [SJ By way of example, Tan's data yield a 
centroid at 2129.55 ± 0.18 MeV and the peak width T = 12.0 ± 0.46 MeV. The 
centroid is higher than in the analysis of Tan who reported 2128.7 ± 0.2 MeV. 
Also the width is much larger than the 7.0 ± 0.6 MeV reported by Tan. The 
reason is that he had fitted two Breit-Wigner functions to also account for a 
higher mass shoulder. We proceed by also fitting two Breit-Wigner distribu- 
tions to all data. Then the peaks become smaller and there is good agreement 
between Tan's result and the one here. The agreement within errors shows the 
adequacy of the present method. 

One feature of the data which has been also ignored so far is a shoulder in 
the peak distribution at higher missing masses. In order to take this shoulder 
into account we have fitted incoherently two Breit-Wigner forms to the data. 
These fits are - on a x 2 criterion - slightly superior to those with only one Breit- 
Wigner. A fit on amplitude level did not indicate an interference. The peak 
region after subtraction of a smooth cross section fraction is shown for all data 
considered in Fig. [TTJ 

Also shown are the fits with two Breit-Wigner functions. They account also 
for the higher mass range which was not in the case with a single Breit-Wigner 
(compare Fig. [3]). The fits with two Breit-Wigner functions are also superior 
to those with only one Breit-Wigner on a x 2 basis. This can be seen from 
the compilation in Table [1] The many values with x 2 /dof<l indicate that the 
authors have been very conservative in estimating the errors. 

Obviously, there is a clear peak present in all data sets independent of the 
reaction. The shoulder visible for larger invariant masses appears also in all 
data sets. However, the latter looks somewhat different for the two reactions. 
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Figure 8: Comparison of spectral shapes in the peak region for the inclusive data 
from HIRES [T2] and TOF [33|38]. The cross sections d 2 a / dmdVt* deduced for 
the HIRES data are shown as shaded area indicating the error band. The data 
da /dm deduced from the TOF II data set divided by 300 are shown as full 
circles with error bars. 



Table 1: The reduced x 2_va hies obtained in the fits wit one Breit-Wigner (1 
BW) and two Breit-Wigner forms (2 BW). 



Ref. 


xVdof (1 BW) 


X^/dof (2 BW) 


Braun 


0.64 


0.63 


Eastwood 


0.92 


0.58 


Pigot 


2.93 


1.62 


Tan 


0.97 


0.41 


SPES4 


1.5 


1.2 


HIRES 


0.42 


0.32 


TOF II 


1.63 


1.0 



The fit with two Breit-Wigner functions accounts well for the data sets in all 
cases. 

The fit parameters for the two Breit-Wigner forms are visualized in Fig. 
[T2l Again the values obtained for the data set TOF I are larger than all other 
results. The mean values and variances are also shown and given in Table [3J 
The values for the x 2 /dof indicate rather poor fits. For the lower-mass peak 
the centroids moi scatter around the result for Tan's data. In all cases the fit 
parameters scatter. A plausible reason for this could be uncertainties in the 
beam momenta. However, then one would expect both centroids for a given 
measurement to be above or below the mean value. This is not the case. It 
should be mentioned that the second higher mass structure for the pp — > K + Ap 
reaction is smaller than for the K~d — > n~Ap reaction. This may point to 
differences in the reaction mechanism for the two processes, which has to be 
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m(Ap) (GeV) 

Figure 9: Data after subtraction of cross section with smooth behavior (for the 
K~d — > 7r~ Ap reaction) or subtracted phase space distributions for pp — ¥ K + Ap 
and pp K + Yi + n as discussed in previous sections. The data (full dots with 
error bars) are from Braun et al. [5] , Eastwood et al. [5] , Pigot et al. [TU] , Tan [3] , 
SPES 4 [TT], HIRES HI], and TOF [37]. The data from gUnHHHD] are from 
strangeness exchange $2$, the others are from associated strangeness production 
0. Also shown are fits with a Breit-Wigner distributions (solid curves). 
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Figure 10: The fitted centroid values (lower panel) and widthes (upper panel) 
from the fits with one Breit-Wigner distribution for the data sets from Braun 
et al. ©, Eastwood et al. [9 , Pigot et al. [Tfl], Tan % SPES4 [11], HIRES [12], 
and TOF [37]. The lines show the mean (thick line) and its variance (thin lines). 
The arrows show the positions of the indicated thresholds. 



expected anyway. 

5.1 The lower-mass peak 

We are now left with the question of what the origin of the structure is. We 
have shown above that the lower-mass peak position coincides with the E + n 
threshold. This would hint to a genuine cusp effect |14j . However, the question 
arises to what precision the experimental missing masses are known. Possible 
shifts in the missing mass distributions in the bubble chamber measurements 
were minimized by studying the chamber magnetic field, energy range relation 
or kinematic fits to reference reactions. In the case of the kaon beam applied 
by Pigot et al. [10] the incident beam was carefully measured with the help of a 
magnetic spectrometer. In the proton-proton experiments the question reduces 
to the precision with which the accelerator beam momentum is known. In 
lower energy experiments at COSY Jiilich employing the spectrograph Big Karl 
deviations of 0.11 to 0.5 MeV/c were found [35] and the COSY beam energy 
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Figure 11: Same as Fig. [9] but for two Breit-Wigner distributions fitted to each 
data set. 



in later experiments was corrected for these small deviations. Similarly for 
momenta around 1930 MeV/c a deviation of 2 MeV/c toward higher momenta 
was found [30]. Such a deviation by 2 MeV/c would amount to a shift up 
of 0.7 MeV for the HIRES experiment. However, no correction was found to 
be necessary when taking the Ap threshold as benchmark 35 . Note that the 
demanded agreement between the HIRES and the SPES4 data required a shift 
of the latter to the low energy side closer to the nominal accelerator beam 
momentum. From these findings we conclude that the energy scale in the data 
used here is reliable and therefore the lower-mass peak is at the E + n threshold 
within two standard deviations of the experimental uncertainty. Thus, there is 
no evidence for a resonance in the Ap system below the EiV threshold which in 
turn restricts the possibility of the existence of a strangeness S = — 1 counterpart 
of the deuteron [Hj . Of course, nothing can be said for the region within that 
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Figure 12: Same as Fig. [5] Left panel: The Breit-Wigner parameters for the 
low mass structure for the data sets from Braun et al. [Q , Eastwood et al. [9], 
Pigot et al. [ID], Tan 0], SPES4 [IT], HIRES [T2], and TOF [37]. The lines 
show the mean (thick line) and its variance (thin lines). The arrows show the 
positions of the indicated thresholds. Right panel; same as left panel but for 
the high mass structure. Note the wider scale on the right side compared to the 
left side. 



uncertainty so that a quasi-bound deuteron-like state extremely close to the 
threshold cannot be excluded at present. We will study this possibility 
further. 

5.2 Flatte approach 

At this stage let us emphasize that the Breit-Wigner distributions are used here 
solely as a tool to determine the peak positions of the structures seen in the 
various experiments and to obtain values for the widths that can be compared 
with each other. They are not meant as a physical interpretation of the data. 
There are other parameterizations that are considered to be more adequate for 
amplitudes near a threshold like the one proposed by Flatte [12] , specifically, if 
one wants to determine also the pole positions. However, those forms usually 
involve also more free parameters. Furthermore, the finite resolution of the 
experiment requires that the yield should be folded with Gaussians so that the 
functional form in the peak region would be identical anyway (see Braun et 
al. i). 

In any case, we performed also exploratory calculations with the relativistic 
formulae of Flatte [12] • That the peak in Ap occurs just at the T, + n threshold 
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Table 2: The mean of the fitted Breit-Wigner distribution parameters excluding 
the TOF I data set. The centroids are denoted with moi and the widths with 
IV 



variable 


fitted value (MeV) 


xVdof 


m i 


2128.7±0.3 


8 


Ti 


8.0±0.7 


6 


m 02 


2140.0±1.4 


23 


r 2 


11.3± 2.0 


6 



is no accident, if the production of the two hyperon-nucleon states are viewed 
as occurring via s-wave hyperon-nucleon scattering. As Flatte has pointed out 
then the drop on the high-mass side may be interpreted as the effect of the 
opening of the S + n channel operating through unitarity. But the imposition 
of analyticity requires that the presence of the S + n channel must be felt below 
threshold as well, thus creating the rapid decrease in the Ap cross section on 
the low-mass side of the peak. The differential cross section around the S + n 
threshold is then 



da v 
dm 



C 



m 2 — im r (Ti + T 2 ) 



(7) 



with v = 1 denoting Ap Ap and v = 2 denoting S + n — > Ap. The partial 
widths are 

I\ =T{Ap^ Ap) = g 1 k 1 (8) 
with gi an (effective) coupling constant squared |42j and k\ the Ap cm. momen- 
tum so that m = m\ + k\ + \Jrri 2 , + k\. 



c?. C denotes an overall normalization 



constant which includes the production cross section. m r is the mass of the 
resonance and Tq the partial width into the incoming channel. It is usually 
assumed to be weakly momentum dependent and thus can be assumed to be 
constant 42], For the other channel we have 



T 2 = r(E+n -> Ap) = g 2 k 2 



(9) 



with fc 2 the S + n cm. momentum. Above the S + rt threshold k 2 is real but 
below the threshold it becomes purely imaginary so that 



T 2 = r(S + n -t Ap) = g 2 K 2 
with n 2 — i\k 2 \. In this region the cross section is given by 

m^r ri 



rfcr i = c 

dm (m 2 . 



m r r 2 y 



(10) 



(ii) 



while da 2 /dm = 0. 

In the actual calculations we assumed that Tq = giki(m r ). In practice, 
To can be absorbed into the overall normalization factor C. It turned out 
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that the fits are not sensitive to the resonance parameters, something that was 
also seen already by Braun et al. [6] when analyzing their own data. Thus, it 
seems difficult to learn more on the pole position. Therefore, we kept those 
parameters fixed to reproduce a peak at the (S + n) threshold and fitted only 
the other parameters (couplings and an overall normalization). In addition we 
assumed a Gaussian smearing with a fitted width to take into account the finite 
resolution of the different experiments. Although there is a large variation in the 
obtained values, the ratios of the coupling constants arc better defined. Indeed, 
this is just a reflection of the scaling properties of the Flatte parametrization 
discussed in Ref . [55] . It is worthwhile to mention that the final curves due to the 
folding with the Gaussian resolution function arc quite similar to Breit-Wigncr 
functions. 

The standard Flatte parametrization Eq. ((7J| produces only a single peak. It 
cannot generate the secondary peak (or shoulder) seen in the data. In order to 
account for this component we have coherently added a Breit-Wigner function. 
However, in the actual fit procedure the centroid of the Breit-Wigner is always 
moved to the maximum while simultaneously the Flatte part is reduced. The 
fits now show a wider distribution but with almost no interference. It should 
be stressed that the assumed Flatte formalism is only for a coupled two-channel 
system. In the discussed AN— UN system there is also a tensor coupling between 
the 3 5*1 and 3 Di partial waves so that effectively, one deals with a system of 
four coupled channels and one would need to use a suitably generalized form for 
the Flatte parameterizations [3]. 

5.3 The higher- mass structure 

In the case of the second structure, which in most of the data sets looks like 
a shoulder, there is no near-by threshold. Therefore one possible explanation 
is indeed that there could be a genuine resonance in the Ap system. However, 
such structures could be also artefacts of the cuts applied to the data from the 
K~d — » 7r~ Ap reaction in order to exclude events from processes with the proton 
being a spectator. The momentum distribution of the nucleons in the deuteron is 
asymmetric with respect to its mean, say a Hulthen function with a maximum 
at 50 MeV/c. Energy cuts on the proton spectra in this order of magnitude 
were applied as discussed in Sect. [3] Deloff [13] showed that such cuts on the 
proton energy as well as on the proton emission angle affects the spectrum for 
m (Ap) < 2.1 GeV. In order to study whether such cuts can produce structures in 
the spectra at higher mass values we perform Monte Carlo (MC) calculations. 
The quasi-free production is ignored. In addition to phase space the lower- 
mass peak is represented by a Breit-Wigner distribution ([6]) with too = OToi 
and r = Ti the corresponding parameters from Table [5] A beam momentum 
of 760 MeV/c was chosen, which corresponds to the case of Ref. [6]. The Ap 
spectrum obtained in this way is shown in Fig. [13] as black curve. In a next 
step the requirement p p > 150 MeV/c was chosen (red curve). The effect on 
the spectrum is minor. We then introduced the requirement cos(AT~, n~) > 0.9. 
This reduces the height of the spectrum severely with, however, no change in 
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Figure 13: Monte Carlo calculations for the strangeness exchange reaction 
K~d — > ir~Ap at a beam momentum of 0.760 GeV/c. A peak with Breit- 
Wigner shape for the peak with mass mol and width Ti from Table [2] is added 
to phase space distribution. The spectrum without cuts is shown in black, with 
a cut on the proton momentum is shown in red and the one with an addition 
cut on the pion emission angle is shown in blue. 



the spectral shape (blue curve). This is not surprising since below the peak 
there is phase space distribution and hence isotropic emission. 

Another effect studied to inspect variations of the spectral shape are crossed- 
channcl resonances. The Dalitz plot shown by Braun et al. 6 showed the 
excitation of the £(1385) clearly. This resonance at the beam momentum of 
w 760 MeV/c does not have an overlap with the peak at « 2129 MeV. This is 
also visible in the Ap spectrum shown in Fig. 1141 Resonance mass and width 
for that simulation are taken from the PDG [44] . A putative £(1480) resonance 
(referred to as £(1480) bumps by the PDG [44]) is not visible in the projection 
of the Dalitz plot on the ir~A axis in [6]. Also such a resonance would not lie 
below the peak at 2129 MeV. The PDG 44 suggests in a note that there are 
two £ resonances around 1670 MeV: a P n at « 1660 MeV with V « 100 MeV 
with three stars and a D13 at « 1670 MeV with width T « 60 MeV with four 
stars. In order to account for these resonances we have assumed a mass of 1660 
MeV and a width of 100 MeV. Although the Dalitz plot in [6] showed that the 
strengths of these resonances, if excited, are smaller than for the £(1385), we 
assumed the same production cross section, since we are only interested in a 
possible change of the spectral shape. The peak area is clearly shifted up by 
the introduction of a such resonance, cf. Fig. Qj] but the cross section below 
the peak remains smooth. Finally we examine the influence of crossed-channel 
resonances on data taken at kaon beam momenta almost at rest. This is shown 
in Fig. HU 

A possible explanation for a shoulder at higher-mass was given in Toker et 
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Figure 14: MC calculation for the Ap spectrum for the reaction K~d — » tt~ Ap 
for a kaon momentum of 0.76 GeV/c. Influence of crossed-channel S resonances. 
Same description of the calculation as in Fig. [T3J 



al. [23 . In that work a shoulder was observed due to an interference between the 
"direct " A production mechanism and those involving the 'EN — > AN transitions 
amplitude, in case of their interaction that produces an inelastic virtual state. 
We should add, however, that such a scenario was later on disputed by Deloff 
[43] . Rather, he advocated an interpretation that resorted to the presence of a 
P-wave resonance at 2140 MeV. In fact, enhanced P-wave contributions around 
the UN threshold are certainly possible. For example, an inspection of the 
Ap — > Ap partial waves of the Jiilich YN model, cf. Fig. 13 in [15], reveals the 
^i-wave shows indeed some structure. As a consequence, there is a shoulder in 
the elastic scattering cross section at a few MeV above the UN threshold, see 
the solid line in Fig. [2j Its maximum is at 2139 MeV, when transformed to the 
Ap invariant mass, which is in good agreement with the experimental findings 
(see Table [2]). However, the width is certainly larger than the experimental 
results. 

The differences between the data for the two reactions are that the data from 
reaction K~d —> tt~ Ap are total cross sections whereas those for the pp — > K + Ap 
reaction are differential cross sections. Hence in the latter case interferences are 
possible which are absent in the first case. We study the second peak in the 
case of the second reaction further. Here we concentrate the discussion on 
the HIRES case. We may assume the second peak structure as a statistical 
fluctuation. A fit with only one Brcit-Wigner yields x 2 /dof= 22.4/56, while a 
fit with two Breit-Wigner functions yields 15/53. At this point it should be 
mentioned that the large error bars are not due to poor statistics but due to 
systematic uncertainties like flight paths in the spectrometer and decays along 
this path as well as uncertainties in the acceptance. Therefore a good fit does 
not yield x 2 /dof= 1 as is the case of statistically distributed data. In a further 
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Figure 15: Same as Fig. [14] but for only £(1660). The kaon beam momenta 
are next to the appropriate curves. The upper curves (black) are without cuts 
while the lower ones (red) are with cuts. 



step we fit one Breit-Wigner by omitting the data in the range of the second 
peak. The remaining cross section was then fitted by a second Breit-Wigner. 
This is shown in Fig. [T5] together with a fitted Breit-Wigner distribution. Note 
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Figure 16: The remaining cross section after subtracting the low mass peak 
from the HIRES data. The solid curve is a fit with a Breit-Wigner function. 

that the error bars contain the uncertainties from the different subtraction steps. 
The peak shows a significance of 4. Oct. This is smaller than the common practice 
of requiring 5<r for a discovery. 
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6 Summary and Discussion 



We have studied the Ap interaction in the vicinity of the E + n and E°p thresh- 
olds. Experimental elastic transition cross sections do not show any particular 
enhancement at the threshold for the Ap — > EiV scattering. However, the qual- 
ity of the existing data is poor and does not allow for reliable conclusions. While 
the enhancement predicted by the recent Jiilich OBE model [TS] is modest, it is 
much larger in the Nijmegen soft-core potentials [T7HT8] and in the leading-order 
YN interaction obtained within the chiral effective field approach |19) . 

In the study of the Ap FSI in case of the strangeness exchange reaction 
K~d — > TT~Ap the experiments made cuts of 75-150 McV/c on the proton 
momenta, thus excluding real spectator protons. After applying the cuts the 
peaks at the EiV thresholds from different experiments became comparable with 
respect to their shape. The influence of the momentum cuts as well as angular 
cuts were studied in detail by Deloff [43] . 

The enhancement is also visible in high resolution experiments employing 
the reaction pp — > K + Ap. A good agreement is found with respect to the 
shape for those data from strangeness production Refs. [TTjrTJ] and [37]. When 
comparing the latter, where the data from Refs. [11] and [12] are transformed 
into the cm system, with those from the reaction K~d — > 7r~Ap, agreement 
is found for the structure around the EiV threshold with cuts applied. This is 
expected since in the pp — > K + Ap reaction both nucleons participate due to 
7r or K exchange. However, it seems mandatory to augment the pp — > K + Ap 
data when new exclusive measurements with high resolution close to threshold 
become available. 

The coincidence of the peak maximum with the E + n threshold within the 
experimental errors suggests that the peak is most likely due to a cusp effect 
induced by the strong coupling between the AN and EiV channels. Some of the 
Nijmegen YN potentials produce a cusp in the 3 Si — 3 D\ partial wave |18ll2"5] . 
The same is true for the recent Jiilich model [15) . Conclusive evidence for 
an unstable (deuteron-like) bound state in the UN (isospin 1/2) channel, as 
predicted by other YN potentials in the literature [MJ17| , would require a peak 
that is well separated from the (and below the) E + n threshold. 

In both reactions there is a second peak at invariant masses well above the 
E°p threshold. A possible interpretation due to deformation of the lower-mass 
peak by experimental cuts, crossed-channel resonances or interference effects 
can be excluded. Some YN potentials like the recent Jiilich model [TS] predict 
an enhanced P-wave contribution around the EiV threshold, cf. the elastic 1 P X 
phase shift in Fig. 13 of that reference. This leads to a shoulder in the Ap cross 
section in the relevant energy region, which at least qualitatively resembles 
the shoulder seen in the reactions discussed in the present study. But more 
exotic interpretations in form of a genuine Ap resonance are not excluded. For 
example, it might be a S = — 1 state as predicted by Aerts and Dover [46|l47] . 
However, their predicted 1 P\ state lies below the EiV threshold. A recent review 
of S — — 1 dibaryonic states can be found in Ref. [JS] . 

Model calculations that study the enhancement seen in the Ap invariant 
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mass near the SiV are all based on YN interactions that include the coupling 
between the AN and EN channels for obvious reasons [23l|4TJ|49l|50]. The 
investigations by Toker et al. [23] and Torres et al. [31] deal with the reaction 
K~d —> ir~Ap and are performed in the Faddeev formalism. In Ref. [23J different 
potential models are employed that produce an inelastic virtual state (and thus 
a threshold cusp) but also an unstable bound state, where in the latter case 
a resonance-like behavior in the Ap cross section is generated slightly below 
and even above the threshold. No direct comparison of their results with the 
Ap spectrum of Tan [3] is presented. Torres et al. [JT] achieve quantitative 
agreement with the data. They claim that a virtual state near the UN threshold 
is required for reproducing the data but not an unstable bound state. In the 
calculations for the associated strangeness production reaction ((4]) Deloff [49] 
as well as Laget [50] use YN interactions that produce a virtual state, namely 
the models presented in Ref. [33J in the former case and the Nijmegen model 
ND [21] in the latter. 

New experiments with better statistics are required in order to discriminate 
between the different scenarios. A hyperons at low momenta could be produced 
with high intensity at the J-PARC facility in Japan and then used for pertinent 
experiments. Also the strangeness exchange reaction could be studied with the 
E31 setup again at J-PARC [5J. With regard to the pp — > K + Ap reaction 
exclusive data with full acceptance but even higher resolution than the existing 
ones would be highly desirable. 
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